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The ability to detect the interaction of light 
and matter at the single-particle level is becoming 
increasingly important for many areas of science 
and technology. The absorption or emission of a 
photon on a narrow transition of a trapped ion 
can be detected with near unit probability [Tl [2], 
thereby enabling ultra-precise ion clocks [3, 4 J and 
quantum information processing applications |5j. 
Extending this sensitivity to broad transitions is 
challenging due to the difficulty of detecting the 
rapid photon scattering events in this case. Here 
we demonstrate a technique that enables the scat- 
tering of a single photon on a broad optical transi- 
tion to be detected with high probability. Our ap- 
proach is to use an entangled state to amplify the 
tiny momentum kick an ion receives upon scat- 
tering a photon. The method should find appli- 
cations in spectroscopy and quantum information 
processing. 

An ion in an electric trap provides an excellent system 
for carrying out precision spectroscopy. Laser cooling [5] , 
either directly or sympathetically via an auxiliary ion [7] , 
minimises thermal line broadening. Long storage times 
allow for many repeated measurements on the same par- 
ticle and collisional shifts are non-existent. A key require- 
ment is to detect whether photons of a given frequency 
are scattered. For some narrow ionic transitions, cor- 
responding to long excited-state lifetimes, the electronic 
configuration change associated with photon absorption 
or emission can be detected with near 100% efficiency us- 
ing the electron shelving technique [T| . In the general case 
where electron shelving cannot be implemented directly 
on a 'spectroscopy' ion of interest, efficient detection is 
possible using a co-trapped auxiliary 'logic' ion following 
the technique of quantum logic spectroscopy [SJ. Here, 
information about the electronic state of the spectroscopy 
ion is mapped via a joint vibrational mode to the elec- 
tronic state of the logic ion, where it can be read out. 
For broad transitions, quantum logic spectroscopy fails 
due to the extremely short excited state lifetimes and 
resulting inability to spectrally resolve their vibrational 
sidebands [2J. 

Another signature of a photon scattering event is the 
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recoil kick an ion receives upon absorbing or emitting a 
photon [9 1. The size of the recoil can be characterised 
by the ratio of the recoil energy E roc to the energy of a 
quantum of motion hv of an harmonically trapped ion 
oscillating at frequency v in the form of the Lamb-Dickc 
factor r]= ^ E Tec / (hiy) , which typically fulfils 77<gCl in ex- 
periments on optical transitions. For an ion cooled to 
its motional ground state, the probability of being pro- 
moted to the next higher vibrationally excited state by 
the recoil of an absorbed photon is rj 2 . Consequently, 
the method of monitoring the ground-state population 
is inefficient for detecting absorption events. Spectral 
lines from broad transitions have been reconstructed by 
observing changes in the fluorescence of a laser-cooled 
control ion via scattering of many 100's of photons [10] 
or via an off-resonantly induced optical dipole force [llj . 

We now describe a technique for amplifying the recoil 
signal of single photons, which is independent of the par- 
ticular ion species employed (Fig. nl . As with quantum 
logic spectroscopy, a co-trapped logic ion encodes a two- 
level qubit with ground and excited states \\) z and |t)z> 
respectively. The logic qubit and a common vibrational 
mode of the two-ion crystal are prepared in their respec- 
tive ground states. A state-dependent force applied to 
this mode splits the vibrational state into two parts, 
each correlated with a different eigenstate in the logic 
qubit |12[ 113] . Consequently, logic qubit and motional 
state become entangled into a Schrodinger cat state of 
the form $>=(\^) x \+ a ) + \ 4r-) x \-a))/V2, where |->) x 
and \<—)x are eigenstates of the a x Pauli operator and 
| ± a) are coherently displaced motional states. 

In the next step, spectroscopy light is sent in at a 
known frequency. The inverse cat generation opera- 
tion is then applied, thereby recombining the two vi- 
brational components of the cat state and disentangling 
them from the logic qubit. If no spectroscopy photon 
was absorbed, the initial state of the logic qubit | \) z 
is recovered. Otherwise, the recoil of the spectroscopy 
ion causes both vibrational components to be displaced 
by ?7abs so that after recombination the total path en- 
closes an area in phase space. The result is a geonretric 
phase [H] </>abs=2a?7 a bs sin(0 tr ) proportional to the cat 
state size a, which leaves the final logic qubit state ro- 
tated by exp^absCa;) p~5] [16]. Here, (pt I =2irvT where 
r is the time delay between photon absorption and the 
time of largest spatial extent of ^f oscillating at fre- 
quency v. Measurements of the logic qubit spin projec- 
tions (a z ) =— cos (</> a bs) or (a y ) — sm(<p ll \ }s ), using standard 
electron shelving for example, provide information about 



the photon absorption process. In summary, the small 
recoil displacement of the spectroscopy ion is converted 
into a large geometric phase and mapped onto the logic 
qubit's electronic state. 

A complete photon scattering event involves two mo- 
mentum kicks - one in absorption and one in emission. 
Since the absorbed photon comes from a directed laser 
beam, ideally this beam points parallel to the vibrational 
mode supporting the cat state. The recoil direction due 
to the spontaneously emitted photon is random. The 
effect is to add a second, random geometric phase em 
to the measured signal: (<J Z )~ — cos((/) a b s + 0cm), where 
the bar indicates averaging with respect to 4> em . For an 
isotropic emission pattern, this expression reduces to 

(a z ) = - cos Oabs) sine (0 cm ) 

= -cos(2a77 abs sin0tr)sinc(2a77cmsin0t r ), (1) 

and similarly (a y }=sm(<t> abs ) smc(0 e m)- 

In our demonstration we confine a mixed-species two- 
ion crystal in a linear Paul trap (see Methods). The 
isotopes 40 Ca + and 44 Ca + are employed as logic and 
spectroscopy ion, respectively (Fig. ^k). The 40 Ca + 
provides sympathetic Doppler and vibrational ground 
state cooling of the lowest frequency axial vibrational 
mode at v = 1.199 MHz, a long-lived electronic transi- 
tion in which to encode a qubit (S1/2 — > D5/2) an d qubit 
state readout via electron shelving. Spectroscopy is per- 
formed on the strong D 3 / 2 — > P1/2 open dipole transition 
in 44 Ca + , which has an excited state lifetime of 7.1 ns [17] 
and, with 93.5% probability, will absorb only one pho- 
ton before spontaneously decaying to the S1/2 ground 
state PI] , 

The experimental sequence is presented in Fig. [2]d. A 
state-dependent force [Tj|] is realised using a bichromatic 
laser field resonant with the red (—v) and blue (+v) vi- 
brational sidebands of the logic qubit. Precise timing of 
the absorbed photon with respect to the phase of \P os- 
cillating at v is achieved using a phase-coherent train of 
spectroscopy pulses, each of ^60 ns duration and sepa- 
rated by 1/v. Scanning the time delay r between the cat- 
state oscillation and the pulse train gives rise to interfer- 
ence fringes (Fig. 3a). Since the signal (<j z ) (blue curve) 
is insensitive to the sign of the geometric phase shift, the 
fringe period is equal to half the motional period. Con- 
versely, the signal (<j y ) (red curve) is sign sensitive and 
therefore oscillates around with a period equal to the 
motional period 1/v. Furthermore, while the signal in 
(<j z ) is sensitive to the recoils of absorbed and reemit- 
ted photon, the signal in (<j y ) is predominantly due to 
the recoil of the absorbed photon (the reemitted photon 
causes just a slight loss of contrast). The signal in (a y ) 
is therefore, with high probability, due to a single photon 
absorption event. 

A line profile of the single-photon transition in 44 Ca + 
was recovered using the cat state technique (Fig. ^p). 
The ((Ty)-signal amplitude was measured as a function of 
the spectroscopy laser frequency at a reduced light in- 
tensity. The observed line-width of 38(5) MHz is close to 



the width expected due to the combination of the tran- 
sition's natural linewidth (22.4 MHz) with the Zeeman 
splitting caused by an applied magnetic field of 4 Gauss. 
We wish to emphasise that this profile is due, with high 
probability, to only a single infrared photon scattering on 
an open dipole transition. 

The uncertainty in the spectroscopic signal is deter- 
mined by quantum projection noise [20) in estimating the 
expectation values of the logic qubit state. The signal-to- 
noise ratio (SNR) can always be increased by taking more 
measurements. To compare different methods it is there- 
fore useful to renormalise by the number of measurements 
N made and consider the sensitivity (3=SNR/\/N pro- 
vided by each method. For comparison we implemented 
two additional direct detection schemes that do not ex- 
ploit a non-classical state to amplify the scattering signal. 
The first and most basic of the two schemes directly mea- 
sures the probability of the motional state being excited 
from |0) to |1) (phonons) by a single photon scattering 
event by applying a red-sideband n pulse to the logic 
ion followed by a measurement of (er z ). In the second 
scheme, (cr y ) is measured instead and phase sensitive de- 
tection of the motional state is enabled by keeping the 
pulse envelope of the spectroscopy laser phase-locked to 
the beat signal (at frequency v) of the two laser fields 
manipulating the ions' motion. 

Table IT] summarises the experimentally determined sen- 
sitivities. Compared to the direct detection method, 
cat-state spectroscopy enabled an 18-fold improvement 
in measurement sensitivity. Also presented are results 
showing that unlike the direct detection methods, cat- 
state spectroscopy even works without ground state cool- 
ing. Theoretically, the method should work as well with 
thermal states that lie within the Lamb-Dicke regime, 
which can be reached with standard Doppler cooling |21[ . 
Experimentally, a loss in contrast is seen whose origin re- 
quires additional study. 

In principle, the technique's sensitivity can be further 
enhanced by making larger cat states. However, exper- 
imental errors will limit the useful cat size and thereby 
the obtainable sensitivity. In our experiment the domi- 
nant error source is electric field noise on the trap elec- 
trodes, which adds another random geometric phase 0h- 
As a result, the contrast of the signal given in equa- 



tion (HI) needs to be multiplied by a factor exp(— 0^/2) 
with 0h =8 -^h^cat(§T cat +r wait ) where R h is the mode's 
heating rate, n cat = \a\ 2 the cat state size in phonons, 
r cat the time needed for creating the cat state and r wa i t 
the time between cat creation and recombination. The 
larger the cat state, the more detrimental the noise be- 
comes. For the (a y ) measurement in Fig. km and the 
experimental parameters n cat = 8.3(1.0), r cat = 50(2) /is, 
Twait = 32 /is, i?h = 40(20) s _1 , we expect a contrast re- 
duction factor of 0.91(5) and an overall signal amplitude 
of A y — 0.54(4), which is close to the observed value. 

We have demonstrated a method for detecting light 
scattering from broad ionic transitions at the single- 
photon level. Complementing quantum logic spec- 



troscopy for narrow transitions, cat state spectroscopy 
should provide access to a new range of atomic and 
molecular ion transitions for precision spectroscopy. The 
achieved sensitivity could be boosted even further by us- 
ing a larger ion trap to provide a reduced electric field 
noise, or better electronic filtering. Assuming that the 
absorbed and emitted photons have the same wavelength, 
equation ([I]) gives the maximum probability of detect- 
ing a scattering event, in a single measurement, as 61%. 
The limit, imposed by the random emission direction 
of the emitted photon, could be overcome by modifying 
the emission direction, e.g. by employing additional laser 
pulses to stimulate emission, the introduction of a cav- 
ity or by using phased arrays of multiple spectroscopy 
ions. A deterministic measurement of whether an ion 
has scattered a photon would be of interest to the field 
of quantum information processing, especially since this 
measurement could be non-destructive with respect to 
the joint ion-photon state. This desired aspect is ad- 
dressed by recent theoretical proposal by Mur-Petit et 
al. [221. 



Methods 

Experimental setup. Experiments employ a linear 
Paul trap with a minimal ion-electrode distance of 
565 pm. Ions are loaded by isotope-selective photoion- 
ization from a resistively heated oven. The bichromatic 
laser beam used to create and reinterfere the cat state 
is generated by applying two radio frequencies offset 
by 2v to a single acousto-optic modulator, effectively 
producing a beat-note to which the generator producing 
the pulse-train of spectroscopy light is locked. 

Geometric phases contributing to the spectro- 
scopic signal. The total geometric phase that rotates 



the logic qubit state is given by 4> = 4>ahs + 4>cm + 4>h result- 
ing from photon absorption and emission recoil, and ran- 
dom electric fields displacing the ions, respectively. The 
magnitude of the phase (/> a bs=2a?7absSin</>t r is the same 
in every experimental run and is determined by the cat 
state size a and the timing of the photon scattering event 
described by </> tr . The other two phases are independent 
random variables with vanishing odd-order moments. 
This property simplifies the calculation of expectation 
values at t he end of the spectros copy protocol, leading to 

{a z )= - COs(</> abs + 0cm + <M = COsQ a b s )cOS <f > e m COS tjfrfa 

and (a y )=sin(<f> ahs + </> cm + (j> h )= sin(c?!> abs )cos (j) em cos0 h . 
The recoil in spontaneous emission creates the phase 
(f> em =2arj em sin 4>% T cos 9 = (fr cm cos9, where 9 describes 
the angle between the directions of photon emis- 
sion and normal mode oscillation. As spontaneous 
emission from the Px/2 to S'1/2 state is isotropic, 
cos </> em =| /„ d9 sin 9 cos((j>cm cos 9). Carrying out the 
integration results in the sine-function appearing in 
equation ([I]). 

Electric field noise at the trap electrodes with fre- 
quencies close to the vibrational mode supporting the 
cat state leads to displacements of the motional state in 
random directions. White noise acting on a cat of size 
n ca t for a duration r adds a random geometric phase eV 
with zero mean and a variance given by </> b =8i?h"-cat'?" 
|15) . In our model, we assume noise acting on the state 
throughout the experimental steps of cat state creation, 
photon scattering and cat recombination, the steps last- 
ing times r cat , T wait and r cat , respectively. If a spin- 
dependent force of constant amplitude is used for cat 
creation and recombination, we have to make the replace- 
ment r — » |2r cat -I- T wait where the factor | accounts for 
the smaller average cat size during the first and last steps 
(see Supplementary Information). 
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FIG. 1: Cat state spectroscopy in phase space, a) Schrodinger cat state preparation. A qubit encoded in a logic ion, 
initialized in | \) z (inset Bloch sphere), becomes entangled with a joint vibrational mode of the two-ion crystal formed with the 
co-trapped spectroscopy ion. b) Absorption of a photon by the spectroscopy ion causes a displacement of size r; a b s (magnified 
for clarity) in a direction determined by the event's timing relative to the cat-state's oscillation, c) The cat-state is reinterfered, 
disentangling the internal state from the motion and leaving the geometric phase 4>a,bs in the logic ion's internal state, where it 
can be read out via standard electron shelving. 
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FIG. 2: Experimental details, a) Relevant electronic energy levels, b) Experimental laser pulse sequence. The 40 Ca + 
dipole transition at 397 nm is used for laser cooling, optical pumping and state detection of the optical qubit implemented 
on its 729 nm narrow quadrupole transition, which is also used for ground state cooling of the joint vibrational mode. The 
* indicates isotope shifts [53] that allow off-resonant scattering from the respective other isotope to be safely ignored at the 
light intensities employed. Dedicated laser beams resonant with the 44 Ca + transitions are used for initialization (397* nm) and 
spectroscopy (866* nm pulse train). 
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FIG. 3: Cat state spectroscopy results. Single photon scattering signals on the broad D 



3/2 — > P1/2 transition in 44 Ca + . 

a) Interferometeric fringes observed by varying the phase difference <j>t r between the spectroscopy pulse train and the cat state 
oscillation phase (see text). Weighted sinusoidal fits (solid lines) yield a peak-to-peak signal amplitude A in the two detection 
bases (<j z ) (blue curve) and (a y ) (red curve). To ensure maximum contrast for the cat state size used, the data is taken at 
light intensities that ensure a single photon scattering probability of p — 1. b) Line profile of the transition taken in the (a y ) 
detection basis. At each spectroscopy laser frequency the signal difference between the oscillation's minimum and maximum is 
measured. To avoid power broadening the light level was reduced significantly, yielding a smaller signal due to p <C 1. 



TABLE I: Comparison of spectroscopic techniques. The sensitivity of various cat state spectroscopy (CSS) techniques, 
and those which do not employ a non-classical state, are compared (see text). Also given are the number of measurements 
required to reach a confidence of three standard deviations (a) in detection of a scattering event. The ** indicate the result 
for a CSS measurement without ground state cooling. The errors are calculated from quantum projection noise. 
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2.7(1.9) xl0 b 


Phase-sensitive (a y ) direct detection 


0.107(10) 


7.8(1. 5)xl0 2 
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0.338(16) 
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0.109(14) 


7.6(2.0) xlO 2 
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Random geometric phases by electric field noise 
causing motional heating 

Motional heating during the cat state generation, pho- 
ton scattering and cat state recombination gives rise to 
additional geometric phases that reduce the contrast of 
the signal. In this online material we show that due to 
heating, a random phase <&h is produced during cat state 
creation and recombination which is given by 

($2) = X ± ncatRhT ( 2 ) 



where n cat = \a\ 2 is the size of the cat state in terms 
of its average phonon number, Rh the motional heating 
rate, and r the length of the two time intervals in which 
the cat state is created and and subsequently put back 
again into the initial state. 

To model the heating, we divide the time interval 
T = 2t into N small intervals of length dt = T/N and 
follow the time evolution of one of the motional compo- 
nents of the cat state. In each of these intervals, fluctu- 
ating electric fields will displace the motional state by an 
amount dpi — dxi + dpi where the projections of d/3i onto 
the x- and the p-axis are normally distributed indepen- 
dent random variables with variance a 2 — ^Rhdt. This 
can be seen by noting that, after time T, the motional 
state displacement due to heating alone is given by 



N N 



i=l 






dxi + dpi 



and the average vibrational quantum number becomes 



N 



(\(3(T)\ 2 ) = }^((dx 2 ) + (dp 2 )) = R h dtN = R h T 

i 

in agreement with the definition of the heating rate. 

In a first step, we will study the effect of motional 
heating on states that are coherently displaced in phase 
space by some classical well-controlled driving force. In 
the absence of motional heating, we assume that the ini- 
tial state |0) is coherently displaced with a (arbitrary) 
trajectory described by \a(t)). We are interested in the 
case where the (transient) motional excitation by the co- 
herent drive is much bigger than the excitation by the 
random fluctuating force causing motional heating. We 
now look at the effect of a heating event (an undesired 
kick of the motional state) at time ti on the state. Be- 
fore this event, the state was coherently displaced by an 
amount a\ = a(ti) — a(t — 0) where a(t = 0) = 0. 
After the event, the coherent displacement is given by 
a 2 = cn(T) — a(ti). The heating event itself displaces the 
state by an amount d/3j. In the absence of heating, the 
final state at time T is given by 

\a(T)) = D(a(T))\0) = e ~ iIm ^ a ^ D(a 2 )D(a 1 )\0) 



where the phase factor cancels the phase arising from the 
multiplication of the two displacement operations. In the 
presence of heating at time ti , the state is modified to 



e- Um(a '^D(a 2 )D(df3 l )D(a 1 )\0) 
e ii m ( a ld^+df3* a2 ) D ^ a2 + d ^ + Oi j| j 

e iImKdft+d/3*a 2 )| a ( T ) + ^^ 



(3) 



i.e., the state experiences a small additional displacement 
and is multiplied by the heating-induced random phase 



= Im(a* 1 dj3i + dj3* a 2 ) 

= Im(a* 1 dj3i) + Im(d(3* a?) 

— Im(al df3i) — Im(d(3i a?}) 

= Jm((ai - a 2 )* d/3t). 



(4) 



If the coherent displacement a(t) is much bigger than 
the heating-induced displacement, the final state is to a 
good approximation equal to \a(T)) multiplied by the 
product of all heating-induced phase factors occuring in 
the interval [0, T]. Equation (41 can be further simplified 
by splitting d/3i into two random variables dui (dvi) that 
are oriented orthogonal (parallel) to (a\ — a 2 )* and each 
have variances a 2 so that 

<f>i = \ai - a2\duisin6i 

where 6i= ± ir/2 is the angle between the orthogonal 
component of d/3i and (ai — a 2 )*- This term ensures that 
trajectories |a(i)) in opposite directions pick up opposite 
geometric phases: as in the case of the two motional 
components in the generation/recombination of the cat 
state. 

After simple substitution, the heating-induced geomet- 
ric phase acquired by the coherent state a(T) is given by 

N 
i— 1 i 

which reduces to 

* 9 = 25^|a(ti)|du,-sin0 < (5) 



in the case a(Tjv) = a(0) = 0. 

For a cat state that is created and recombined, both 
components acquire heating-induced phases ±$2 of op- 
posite sign. Therefore, the relative phase between the 
two components is $/> = 2$ g . We can calculate its mean- 



squared value: 



N 

(^> = 16^|a(t fc )| 2 (^ 



8R h V |a(i fe )| 2 eft - 8R h / |a(t)| 2 dt 



(6) 



For 

a(t) = (^*/ r rt/ : 

the mean squared phase is given by 

16 



t < T 

T < t < 2t 



(H) 



-RhTTl c 



(7) 



integration of eq. (raj) yields the mean-squared phase 



2r 
($£) = 8R h n cat (— +t wait ) 



(8) 



For a random variable X that has a Gaussian distribution 
with (X) = 0, the equality 



which is the result stated at the beginning in eq. &h. It is 
easy to include an additional waiting time t wa u between 
cat creation and recombination in the calculation. For 



(cos X) = exp( 



(X 2 



y/n cat tJT 

a(t) = { ^n cat 

\J n cat (2T + t wait — t)/T 



t < t holds which enables us to calculate the loss of Ramsey 

t <t < t + t wa u contrast of the cat state spectroscopy from eqs.|7]) or (J8l). 

T -r t wa it fv t 2? ZT -\- t wa it 



